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Using effective field theory (EFT) techniques we calculate the next-to-leading order (NLO) 
spin-orbit contributions to the gravitational potential of inspiralling compact binaries. We 
use the covariant spin supplementarity condition (SSC), and explicitly prove the equivalence 



with previous results by Faye et al. in 'gr-qc/0605139 We also show that the direct 



^~; . application of the Newton- Wigner SSC at the level of the action leads to the correct dynamics 

{^' using a canonical (Dirac) algebra. This paper then completes the calculation of the necessary 

■ ' spin dynamics within the EFT formalism that will be used in a separate paper to compute 

the spin contributions to the energy Hux and phase evolution to NLO. 
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INTRODUCTION 



LIGO 



[l], VIRGO y 







and eventually LISA [3] and Einstein Telescope [j] expect to detect 



Q 



radiation from inspiral 
increasingly important 



ing, binary systems, and building templates for these events has become 
5|, |g] . While for late stages of the inspirals numerical techniques are needed 



[ThSI, for the early stages we may rely on the Post-Newtonian (PN) approximation as an expansion 
in small velocities {v/c) <^ 1 (see 



As argued in the literature 



KjI for a review and further references.) 

-14l|. the expectation that black holes may be close to maximally 
rotating in the binaries [15!] produces a significant impact in the accuracy of gravitational wave 
templates. Therefore, if not for detection, parameter extraction justifies the need for physical 
templates that include spin contributions at higher PN orders beyond the leading effects. In 
addition, PN corrections^including spin are relevant for comparison between analytic results and 
numerical simulations 






In order to produce such templates an important building block is the gravitational potential, V, 
responsible for the dynamics of the bodies in the binary system. The leading order (LO) spin-orbit 



and spin-spin contributions to V have been known in the literature for quite some time 



and spin precession has already been observed in some binary pulsars, e.g. 24]. On the other 



hand, only recent 



spin(l)spin(l) 



il 



le next-to-leading order (NLO) spin-orbit [2514271] . spin(l)spin(2) |28l432l| and 



17 



m, 



3314351] contributions have been computed, albeit with rather dissimilar techniques. 



251, 



261 ] within the 'standard' PN 



The NLO spin-orbit effects (2.5PN) were computed in 
approach [10|, namely by obtaining the metric using Bailey-Israel-Dixqn's spin-dependent stress 



energy tensor 



361 ] and then solving Mathisson-Papapetrou equations |37l ]. On the other hand. 



in 27| the frequency of spin precession was directly obtained in ADM coordinates by a 'suitable 
redefined constant-magnitude spin vector' and subsequently a Hamiltonian was derived that was 
shown to be equivalent to the results in [26|. This is a very useful trick that we borrow in this 
paper, although we will obtain the frequency of precession using a completely different approach. 



The NLO spin(l)spin(2) effects (3PN) were first obtained in [2814301] using a new formalism 



based on the application of Effective Field Theory (EFT) techniques introduced in 



, where it was 



coined Non-Relativistic General Relativity (NRGR) , and extended to include spinning constituents 



in 



391]. An EFT framework^ has shown to be extremely powerful in many realms of classical 



gravitational (and non-gravitational) physics. For example, it has been used to compute the 
NNLO spin-independent contribution to the potential (2PN) 4l||: absorption effects for compact 



objects in binary systems and time dependent backgrounds 



gravitational self-force on extended objects 







43]; the electromagnetic [44] and 



45( 1 ; radiation-reaction effects in t 



le extreme mass ratio 



limit [46]; corrections to thermodynamic quantities for caged black holes 



47 



481 ]: deviations from 



General Relativity at the nonlinear level [49|]; and recently applied to the study of cosmological 
perturbations [5u\. 



The radiation sector of NRGR has been developed further in |5ll . |52| |. The purpose of 
this paper is thus to compute the NLO spin-orbit corrections to the gravitational potential^ 
(within our set of conventions), which toget : 



source multipole moments computed in 



51 



ler with the spin-spin results in [28h31 



331 ] and 



521 ] allow us to obtain all the spin contributions 
in the energy flux and phase evolution to 3PN order. This will be reported in a separate paper^ 531 ]. 



A crucial aspect of the calculations dealing with spin in General Relativity is the choice of 
spin supplementarity conditions (SSCs) [5g, |57] and, in a Lagrangian/Hamiltonian framework, 
the resulting (Dirac) algebra in the reduced phase space. This algebra is expected to be quite 
cumbersome in generic SSCs, already in a Minkowski background [59]. Nevertheless, the results in 



2^, 



301 ] were obtained using the Newton- Wigner (NW) SSC at the level of the action, which was 
shown to be a correct procedure in [29] using standard EFT power-counting techniques. As we 
shall show in this paper, the NW SSC also leads to a canonical structure at linear order in spin, 
at least to NLO. This is consistent with the results in [59] (although formally the comparison only 
applies in the test-particle limit.) 

nn 

Subsequently the potential in [28|, 130] was re-derived using a Routhian approach, introduced 



in [29|], where the (covariant) SSC is only imposed at the late stages of the calculation so 



that we may never depart from a canonical algebra, and full 



NLO spin(l)spin(2) Hamiltonian has also been computed in 32] in ADM coordinates, and the 



consistency was found 



3l|]. The 



equivalence with the results in 



281 ] was shown in [3y]. 



^ See kd| for a review on EFT techniques. 

■^ See [54| for an alternative derivation of the NLO spin-orbit potential. 

^ A comparison with the recent results for the NLO spin-orbit contributions to the radiated power in |55l | will be 
performed in j5^]. 



The basic idea of our approach is the systematic separation of the relevant scales (Afs ^ Ap ^ X,.) 
in the two body problem; from finite size effects (Afs ~ Gm for compact objects), to the potential 
(Ap ~ r) and gravitational radiation (Ar ~ -) scales (for a review see [60(].) For the theory of 
potentials the NRGR prescription consists on computing all possible (Feynman) diagrams (without 
external radiation gravitons) scaling with a definite power of v. For spinning extended objects, the 



31 



33 



3S, 



391]. The EFT 



'rules' are derived from an action (Routhian) principle and are listed in [c 
formalism is therefore significantly more efficient than the (more traditional) FN methods since it 
maps complex integrals into the computation of Feynman diagrams, uses textbook regularization 
techniques, and it is especially suited to handle spin degrees of freedom. Feynman diagrams allow 
for a very natural systematization and (physical) visualization of the computation, and moreover it 



can be automatized almost entirely using Mathematica code |6ll |. As a result, the calculation of the 
NLO spin(l)spin(2) potential was no more involved than obtaining the Einstein-Infeld-Hoffmann 
Lagrangian as in [38], contrary to the intricate approach of j32l ]. 



The NLO spin(l)spin(l) potential (also at 3PN) was computed in 33|. Most of the 



spin(l)spin(l) effects are due to (self-induced)finite size contributions [39] (and an extra term 



responsible for the preservation of the SSC 



29 



33f].) In the spirit of EFTs these terms are encoded 



in higher dimensional operators in the effective action (constrained only by the symmetries of the 
theory) up to unknown coefficients that are 'matched' with a given observable in the complete 



theory 38!], for example the quadrupole moment of a spinning black hole 33|, l39|. Once again 
the gravitational potential follows as a sum of Feynman diagrams of definite FN order which 
greatly simplifies the treatment of finite size corrections. These results were corroborated in 



M, 



m- 



The Hamiltonian of 



34 



35l ] was partially obtained by imposing Poincare invariance, which is not 
sufficient to determine all the coefficients, and the extra terms were fixed 'by a three-dimensional 
covariant ansatz for the source terms of the Hamilton constraint.' As argued in [34] : 'the arguments 
that lead to a unique fixation of the coefficients are, however, far from being straightforward.' 

On the other hand, the determination of the potential in [33,] is straightforward once a (very 
natural) higher dimensional operator is included. Moreover, using EFT power counting techniques 
one can easily show that one (and only one) operator is sufficient to NLO 33, |39|, whose overall 



(Wilson) coefficient is directly related to the spin quadrupole of the Kerr black hole. 



Truth be told, the advantage of the Hamiltonian in 



341 ] is that the spin tensor is fully reduced, 



and one only deals with a three- vector with canonical brackets. Li our formalism we bypass the 



complications of the (Dirac) algebra in the covariant SSC by reducing to a three-vector as the 
final step in the calculation. Therefore, in the intermediate steps we have an extra three-vector 
degree of freedom, S''^^* = S^"^, with non-vanishing brackets (see Eq. (j23p ). that is only solved for 
via the SSC once the equations of motion are obtained (see Eqs. ()22I - I25|) .) This however does not 
represent any major difficulty, no more than dealing with an antisymmetric spin tensor obeying a 
canonical algebra^. 

Ultimately, given the differences in approaches, the agreement found in all cases gives us strong 
confidence the results are correct. In this paper we add yet another piece of evidence for the validity 
of the spin potentials by computing the NLO spin-orbit effects using the EFT framework of NRGR. 
The computation proceeds systematically as a sum of Feynman diagrams. Hence first we sketch^ 
the ingredients that determine the necessary Feynman rules before we embark in the computation of 
the diagrams, that we divide into terms linear and quadratic in G (Newton's constant.) Divergences 
are handled by textbook regularization techniques that result in finite integrals. After the potential 
is obtained we perform a spin redefinition that takes the spin dynamics into precession form from 
which we read off the precession frequency. The details are summarized in an appendix. Then it 
is straightforward to show that our result resembles those in [26|| since both frequencies, ours and 



the one given in 



26t |. are related by a total time derivative. Following [27| one can then construct 



2gl, and 



a Hamiltonian whose equations of motion turn out to be canonically related to those in 
the equivalence is thus proven. 

In the last section we show that imposing the NW SSC at the level of the action also leads to 
the correct dynamics in a canonical manner. Remarkably, the transformation found in [58] between 
covariant and NW SSCs in a flat background applies to our case once re-written in a locally fiat 
frame and with derivatives taken with respect to the 'free falling' proper time. Therefore, here we 
extend the results in \5^ to the case of self-gravitating objects, as long as we ignore spin(l)spin(l) 
(finite size) terms^. 



■* One can also perform a transformation that reduces the spin 3- vector algebra to standard brackets [35| . 

^ We encourage the reader to consult the NRGR literature for more details [2^, [29|, [3l|, [SJ, ISa . |39| . |6(J | . 

^ Imposing the NW SSC does not lead to canonical brackets in the spin(l)spin(l) sector since it knocks off a 

necessary extra term in the Routhian, proportional to the SSC, which guarantees its preservation upon time 

evolution HIs^,!!!. 



II. SPIN EFFECTS IN THE EFFECTIVE FIELD THEORY APPROACH 



2^, 



b include spin effects in the EFT formalism the following action/Routhian was introduced 



29 



31 



33, 



39| 



a=l,2 V ^ / 



(1) 



where the dots include 0(8^) corrections that are not relevant for the present paper^ 
The equations of motion derive from 



IJnax^,, '-^^is'^ni ^M^.s-). 



(2) 



where the potential is given by y = —TZ and 



|^ab^^cd| ^ ^acgbd^^bdgac_^adgbc_^bcgad^ 



(3) 



Following the NRGR prescription |38|, |39| we need to compute all possible diagrams linear in spin 
that enter at 0{v^) in the gravitational potential. We use reparamaterization invariance to select 
X = x^ = t such that v^ = -^ = (1, v). Then each diagram contributes a term in the effective 



action/Routhian which is related to the potential via i J dtTZdia.g 

m, y, y, y 



The necessary Feynman rules are 



: / dtVd 



lag 



m, 



31 



33, 



38|, 



m- 






^Sg 



^sg 



-^sg 



2mr. 



-HiokS 



ik 



{HijhS^ W^ + HnahS^ 



2mr. 



lij,k 



Ok 



— ( HQjhS v-' + HionS^ 



2mr. 



2m, 



kO^i 



-HikflS V*, 



(4) 
(5) 
(6) 
(7) 



A similar expression was advocated in [63] • 

In the approach of (29|, l3ll . 133] we were able to sweep most of the technicalities of the preservation of the SSC into 
an extra term in Eq. ^ proportional to spin(l)spin(l) that we can ignore at linear order in spin and also for 
spin(l)spin(2) effects as it was shown in [281. ISOl. ISlj. 



for the terms linear in the metric, whereas to 0{H'^ 



L" 2 = 



S^^ H, {Hox^i — Hoi^x) 



^^' " Ami L' 



S"^^ Hj {Hki,i - Hkid) V + S"" lHooHoo,i + H^Hqq^i 



(8) 
(9) 



where we only kept the terms that lead to non-zero contractions at the desired order. The field 
Hny represents the potential graviton responsible for the gravitational binding of the binary system 
38|. 

To compute spin-orbit effects we also need the Feynman rules for the spinless couplings. These 
are 13811 












LI 



m 



2mr. 
m 



nir. 



m 



2mp 
m 

2mr. 



-H, 



00 



I0i> 



HtjvW^ + o^oov^ 



i^0ivV2 



(10) 

(11) 

(12) 
(13) 



and for the non-linearities 



TV- 



tV 



m 
8m| 

m 
2m| 



HooHoo 



lOO^Oi' 



(14) 
(15) 



In all these expressions rr ip = , „ ■ We encourage the reader to consult the literature for more 



detailed discussions 



28|, 



m 



31 



Vl387r« 



33|, 



38|, 



3S 



6C| 



A. The (covariant) spin supplementarity condition 

In the Routhian formalism described above the (covariant) SSC is only enforced once the equa- 
tions of motion are obtained using Eq. ([2]). However, at NLO we need to account for corrections 
in the spin algebra due to gravitational effects in the transformation between local and PN frames. 
Since we work with the spin in the local frame, the covariant SSC has the form 



QiO _ aij^l 



(16) 



where 

Gm2 \ 2Gm2 



vr = vl^l + — ^J-— -^v^, (17) 

vr' = (i-^ (18) 



is the particle's velocity in a locally flat frame, v"" = e"u^, and 



eg(xi) = l-^ + --- (19) 

eg(xi) = -2^^ + ... (20) 

e)(.,) = S)(l + ^^+... (21) 



3ii 



with r = xi — X2. The spin algebra in Eq. ([3| then takes the form 

{S\ S^} = -e^^■'=S^ and {S^\ S^} = -e'^'^S^'' (22) 

or equivalently^ 

{S^S(°)^■} = -e^^■^S(°)^ (23) 

where we introduced the spin (three-)vectors S* = e^^ S^ and S'°^* = S"*^*. 

Even though we have to deal with S^^^ all the way to the end, it is straightforward to derive 
the equations of motion, for example for the spin dynamics, using Eq. (I22p . If we write the S and 
S*-"-* contributions to the gravitational potential as 

Vs(o) = A • S(°) , and Fg = ^s • S (24) 

respectively, then using Eq. ([2]) the spin equation of motion becomes 

f =-sxS+(jxs)xA (25) 



^ Notice that our conventions differ by an overall sign with those in [35 



10 



in the covariant SSC, namely 



S(o) 






(26) 



Notice that the IPN 0{G) corrections in ([m [H]) induce 0{G'^) terms at 2.5PN in Eq. ([25]) 
after inserting the LO part of A (that is, say for particle one, A^^ = a^^ = —Gm2T/r^), similarly 



to what happened in the spin-spin sector 



30, 



31 



33| 



III. THE NEXT-TO-LEADING ORDER SPIN-ORBIT POTENTIAL 

To compute the spin-orbit potential we proceed systematically. First we compute diagrams with 
a single graviton exchanged, as in Figs. Hh-/. Then we move to the non-linear gravitational effects 
that can be split in two groups: the seagull diagrams depicted in Figs, ^-d that follow from Eqs. 
(0 [9]) and (J14l llSp . and the 'three-graviton' interactions shown in Figs. [Sli-g, where we contract 
the worldline couplings linear in Hfj_^ with the 0{H^) (bulk) interaction from Einstein's action (in 



background harmonic gauge) [3a]. Everywhere in this paper J stands for J ,^3^ , and j dfs are 



suppressed. The results for the diagrams are collected below. 

A. One-graviton exchange 

It is sometimes convenient to group terms together while computing the one-graviton exchange. 

For instance, we can obtain most contributions in Figs. [TJx-c by (Wick) contracting the 0{H) spin 

li 

coupling [39| 



H^^^pS-Pv^^, 



(27) 



with the spin-independent part 



"" H^pV^vP, 



2m. 



(28) 



and using the two point function 3S ] 



(F^,(x?,xi)i?„^(xO,X2)) = -i5{x\-xl)P^p^, [ ^e^P" 

Jt) P 



(29) 



11 



where 



Pabcd = 2 iVacllbd + Vadribc - VabVcd) ■ 



(30) 



The remaining piece conies from a contraction between the terms in Eq. ([5]) and the one propor- 
tional to Hqqv'^ in Eq. (fT^ . Hence the result for these diagrams is 



m. 



Gmo 



be 



r' \ Sf ( ^^rl - 2vi • V2 ) + 2Sl'viv, • V2 - -^S\'v{ \ + Sf ((vi - 2v2)V2 • r - 2a^r2) 

(31) 



while the acceleration-dependent piece results from the (partial) time derivative (p = 0) in 

Eq. (p7|) . The latter induces derivatives of the particle's position and velocities, as well as its 

spin. In order to avoid contributions that depend on Sf, we evaluate them as explained in 

31 1 (see Eq. (46) in 31(1)) namely by using integration by parts so that we may set Oq Sf^{xi) = 0. 



The diagrams in Figs. [Hie account for 'propagator' corrections to the LO potential. These arise 
after we include the non-instantenous part of the propagator in Eq. (|29l) . and expand in powers of 



Po 



m, 



38| 



pg-p2 



p2 V P^ 



(32) 



Again we face temporal derivatives that we handle as we sketched above (that is the reason we 
generate a2 and a2 terms, see below.) Then we have 



%e 



Gm2 



5f 



1 



+ r'Sl'ai + 



! 2 2r2 
3(V2 • r)2 



(V2 • r)^ 



1 



-a2-r ] + 



_L(v2.r)2-ia2-r)SlM+25Mv2-r 



^iO 



v2 + as • r 5f v^2 r^ + Si^v^v{v2 • r - r^Sl'4^^i - ^S^^r{ac!, 



i„2 



+ Sr v^V2 • r + -a*2^ 



(33) 



Finally there is the diagram in Fig. [T]/ which is straightforward to calculate. 



2Gm2 ^ji I j 



(34) 



12 



Therefore the total for the one-graviton diagrams reads 



V^ 



Grrio 



Igrav 



5f (2vi-2vi.V2-^(v2-r)2-ia2-r] + ('2vi.V2 + ^^^^-2vi+a2-r)5;V^2 



1 



2r2 



V2 • r)^ + -a2 • r 5}V{ + 25}^a^V2 • r + r'Sl^ki } r' + Sf (vi - V2)V2 • r - -a^r 



+5iV^vjv2 • r - r^S'^SL^vl - -r''s{^aiv\ 



J„2 



(35) 



V 

-0- 



"tir 



v 
O 



V" 

a) 

3 




V 

-0- 



V" 

d) 



e) 



V- 
-O- 



f) 



FIG. 1: One-graviton contributions linear in spin at 0{v^). A dashed line represents a potential graviton, a 
blob account for spin insertions at each PN order from Eqs. IM'\T!), send a black dot stands for 'propagator' 
corrections as explained in the text. 



B. Non-linear gravitational interactions 

1. Seagulls 

These diagrams are shown in Figs. [JJi-d and account for the 0{H^) spin and spin-independent 
couphngs in the worldhne action/Routhian. The results are 






2^2 



G^m 



G'^m2Tni 
% = 2 -i 



G 7712^-1 






(36) 
(37) 
(38) 
(39) 



13 



yi ^5 ^3 

-O O 7^ 



v° a) v^ v° b) v° v° c) v'- 



v° d) v^ 



FIG. 2: Non-linear gravitational effects due the couplings in Eqs. HI [3) and i[74l[l5j) . 



ich add up to 




•'seagulls - 


r4 



mi 



5f + (mi - 2m2)S{^^r{ + 2(m2 - mi)S{^w(j r' (40) 



2. Three- graviton couplings 

Finally we have the three-graviton couplings in Figs. [3]i-(7. These are easily handled with the 
aid of Matheniatica code |6l|. The results are: 



m = -^s--'-2 (41) 

m = ^^i'r^-2 (42) 

% = -^(25fr^ + ^5iVvl) (43) 

% = -S^^^^rVl (44) 

% = ^^rrVl (45) 

Hf = ^^ (-25f r^- + Is^r^-i) (46) 

m = -2^^5rrV^, (47) 



14 



V 

-o- 



V 

-o- 



O- 



a) v^ 



i V 



u" c) u" 



Or 



t^^ d) t' 



.,0 .^ ^ 



e) t.'^ 



U° /) V^ 



V° g) ^ 

FIG. 3: Non-linear gravitational effects due to the three- graviton couplings. A black dot represents an extra 
power of V from a time derivative in the three graviton vertex (see text.) 



for a total 



V, 



G^m2 



Sgrav 



-2(mi + m2)5r + 



^^iM + (^ + 2-,)^iV,)r^ 



(48) 



It is easy to foresee that many three-graviton diagrams entail the following integral 



Iijl(.T^) 



p«piq' 



q,pq^p^(p + q)^ 



q,k(p + k)Vk2 



(49) 



where two powers of momenta come from the three-graviton vertex, and an extra one from the 
spin coupling(s). To compute Iiji we simply reduce it to scalar integrals, and using 



1 



1 



qq2(p + q)2 8|p| 



(50) 



we obtain 



I, 



ijl 



647r2r6 



•2(rM^'+r^5^M-2rVrM . 



(51) 



The values for each diagram follow by contracting liji with the spin tensor and velocities. Those 
involving S^^ are indeed proportional to (derivatives of) {HqqHqqHqq), and can be easily computed 



15 



from the results in 



M- 



Aside from the different powers of momenta, the only difference comparing with the kind of 



diagrams we encountered before in 3l|, |33|] for the spin-spin potentials is the appearance of time 
derivatives in the three-graviton vertex, depicted as a black dot in Figs. [3t and[3b. For example, 
we have terms of the sort ~ HoQdiHQidoHQQ from Einstein's action which produces a P2P3 type of 
Feynman rule, etc. To handle these diagrams we proceed the same way we compute propagator 
corrections"^'^, namely p^ ^^ p • v, which induces an extra power of v in the otherwise (naively) 
0{v^) diagrams in Figs. [3) and [3^. 



C. Divergences 



In the computation of the three-graviton contributions we encounter divergences. These can 
be easily understood diagramatically, as in Figs. Wp-b, and arise once a term proportional to p^ 
(a = 2, 3) from the three-graviton Feynman rule cancels out with one of the (two) propagators 
(~ -^) that couples the vertex to the worldline. For example we have Im = in Eq. (|5ip since it 

Pa 

produces purely divergent integrals proportional to /j^ j^ and /j^ j^, as shown in Figs. Hh and [It 



respectively. These divergences are 'self-energy' type of terms that are absorbed into the coup 



of the theory, and are set to zero in dimensional regularization (see the discussions in 
for more details.) 



M, 



i ngs 



38 



39| 



a) b) 

FIG. 4: Type of divergences in the three-graviton diagrams 



^° For the three-graviton diagrams in Figs |3}) and [3^ we can ignore the extra terms that are proportional to time 
derivatives of S'-' , since these are sub-leading (recall S'-' ~ u^S*-'.) 
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D. Assembling the pieces 



Adding up all the ingredients the NLO spin-orbit potential reads 



K 



NLO 



Grw} 



5f (2vi-2vi.V2-;^(v2-r)2-^a2-r) + (2vi.V2 + ^^^^^-2vi+a2-r)5lV^2 



_3_ 

2^' 



(V2 • r)2 + ^a2 • r S'^^^{ + 2S'^ 3.1^2 ■ r + r^Sl^a^ r* + Sf (vi - V2)V2 • r - ^a 



2r2 






J„2 



Smg 



-(mi + 2m2)5f + (mi - ^) SJ^vl + ^-^S^^i 



(52) 



We can now replace the acceleration terms by using the LO equations of motion, namely a2 
—^r, and including also the LO part of the potential 391], we finally obtain 



r 



so G'm2 



5f ( 1 + 2vi - 2vi • V2 - ^(V2 • rf - ^(3mi + 2m2) ) + 



1 - ^(^2 • r)2 + ^ (4mi - m2)) S^^i -(2- 2vi • V2 - ^^^' ' '^ 



G 



,2 +2v^-^(2mi+5m2))5lV^2|>r' 



+ 5f (vi - V2)V2 • r + 5lV{v^V2 • r 



+ 10 2, 



(53) 



which can be written as (see Eq. ([2 



q=l,2 



(54) 



with (for particle 1) 



Ai 



Gm2 
Gm2 



3 G 

-1 - 2v^ + 2vi • V2 + :^(v2 • r)2 + — (3mi + 2m2) ) r - (vi - V2)v2 • r 

/ 3 G 

-(v2 • r)vi X V2 + ( 1 - ^(v2 • r)^ + — (4mi - m2) ) Li 



+ ( -2 + 2vi • V2 + ^^^^^""^ - 2vi + — (2mi + 5m2) ) L2 



(55) 



(56) 



and Li(2) = r x vi(2). 
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IV. SPIN-ORBIT DYNAMICS 



Using the potential in Eq. ([53]) together with Eqs. ^ [22l [26l) . it is straightforward to obtain 
the spin dynamics, that is given by the expression (see Eq. (|25p ) 






-3^=^Si xSi+(^xSi) X Ai. (57) 



A. Precession equation 

To transform the above equation into precession form first we notice that 

1 1 

Ai = ai + VI X w^LO + -vjal + - vi (vi • aj) , (58) 



with (vi = vi/wi) 



d _ 1 Gm2 

ai = — vi = an + 



3 G 

— (v2 • r)2 - V? + 4vi • V2 - 2vj + -(2m2 + 3mi) ) r 



+ (vi-V2)(2vi-r-V2-r) J , (59) 

ai = ^r, and we used the IPN acceleration from the Einstein- Infeld-Hoffmann Lagrangian 



3a] plus Eqs. (fTTl [T8|) together with the LO precession frequency 

-s^o° = 4? + ^v, X aj = ^ (^Li - 2L2) . (60) 

We may construct now a spin vector with conserved norm with the transformation (see appendix 
A) 

expanded to NLO. The precession equation thus becomes (up to 2.5PN)^^ 

^Si = < X Si (62) 



^^ Notice that the manipulations and results for the spin-orbit frequency are remarkably similar to the results in the 
spin-spin sector [3ll.l33l|. 
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where 



^° = ^Si + ^ (l + ^) VI X Ai = w^o + ^NLO ^ ^^^ X a^ + ivi X Ai + i (vi - vi) x aj, (63) 



with Ai = Ai — al. 



Hence the spin precession frequency takes the final form 






Gm2 



Q 1 g (^ \ 1 

- + -vf + V2 - vi • V2 - -r^{'V2 • r)^ + — (mi - 7712) 1 Li - -(v2 • r)vi x V2 
2 4r^ Zr I 



3(v2 • r)2 o G , , , 

+ ( -2 + 2vi • V2 + ^ % ' - 2vi + — (2mi + Sms) ) L2 



(64) 



B. Comparison with other methods 



The NLO spin-orbit contributions to the dynamics of compact binaries was previously obtained 



m 



2a], with the precession frequency given by 



_ Gm2 
■^BBF — ~~jr' 



+ o'^i + V2 - vi • V2 - -rTV'^2 • r) + 7^ (7mi - m2J Li 



2 8 



2r 



(65) 



-(r • V2) - 3(vi • r) j vi X V2 + f-2 + 2vi • V2 + '^^^'^'^^ -2^1 + — (2mi + Qma) ) L2 



Notice that it remarkably resembles our result in Eq. (I64|) . but not quite. The frequencies differ 



by 



Aa; = .;fBF--^° = ^ 



G G 

— 3?72iLi + 3/71-2 — L2 + 3 (r • vi — V2 • r) vi x V2 



(66) 



which is a total time derivative 



d f Gm2 
Au = —3^- I vi X V2 



dt \ r 



(67) 



This determines the equivalence of both results after the spin transformation 



Si -> Si — 3 (vi X V2) X Si. 

r 



(68) 
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Moreover, we can use Eq. (|64p to construct a Hamiltonian, given by 



-"soC^q 



9=1.2 



(69) 



so that Eqs. (|67l [68|) imply that the equations of motion in 
potential in Eq. ([53]) fully agree [27[. 



26l | and the ones that derive from the 



V. THE NEWTON- WIGNER SSC 

Let us consider now the application of the NW SSC at the level of the action. The constraint 
reads 



oiO 
1 ~ ^0 



~0 ~j 

Pi nij Pi 
1 ^' 



p^ + m Pi 



(70) 



Notice this expression is reparameterization invariant, so we have 



q/^ = = h 

2 dX 



(71) 



which we can replace back into the SSC 



5f 



^0 



v^ 



^1 qij Zl 



(72) 



Then choosing once again A = x'^ we obtain (recall S^*^^* = S^^ 



j(0) 
'l(NW) 



1 + Wl 



VI 



2 ,Si(NW)X^ 



(73) 



We can now apply Eq. ([73]) into Eq. ([M]) and obtain (suppressing the 'NW in the spin vector^ 



TASO 



E 

9=1,2 



1 + A 1 



^xA,)+^s, 



>S„ 



(74) 



from which using canonical (Dirac) brackets, that is (say for particle one) w^^^v ~ ~as^' ^^ S^t 
the frequency given by Eq. ([63]) (recall vi = vi/wj'), and ultimately Eq. ([6l]) including the one 
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'extra term' from the 1/(1 + vT— vj I factor in Eq. ()73p given by 



;Vi (vi X a^) , 



(75) 



after expanded to O(f^). 



Let us make a pause and remind ourselves of the form of the transformation between covariant 
and NW SSCs in Minkowski space given by Eq. (B13) of 



58] 



m 



Snw = -r;-S + 



pS 



Ho Ho{Ho + m) 



(76) 



with Hq = y'p^+m?. Using Eq. ()7ip we have (A = t) 

p = m(lH jv+..., and p = m (v + v +...), 



(77) 



that we can insert back into Eq. ()76p and we get 



>NW 



1 



1 



1 - — - -v^ + . . . S + -v(v • S) 1 + — + 



(78) 



which is precisely the transformation in Eq. (j6ip after the replacement v — t- v. On behalf of the 
equivalence principle, it is indeed not surprising that the same redefinition as in Minkowski space 
still applies if rewritten in a locally flat frame and with derivatives taken with respect to proper 
time (as long as we ignore finite size effects.) 



Our result thus provides (yet another, see J29l ]) explicit example where the NW SSC leads to 



standard brackets, this time at linear order on the spin. Hence the Hamiltonian in Eq. (j69p can 
be equivalently obtained by applying the NW SSC from Eq. (j73p directly in the expression for the 
potential in Eq. (fMl) . such that (see Eq. ([7i]l ) 



Hs, 




^s, + 



1 



1 + Jl 



(Vg X Aq) 



•s„ 



(79) 



VI. CONCLUSIONS 



Using the EFT formalism of NRGR 



38( for spinning extended objects [29|, |39|], we computed the 



NLO contributions to the spin-orbit gravitational potential of inspiralling compact binaries. We 
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explicitly proved the equivalence with previous results in the literature [26|, |23] and provided extra 
support for the canonical aspects of the NW SSC (that was already shown to be the case for the 
spin(l)spin(2) sector to 3PN order in [29].) This is consistent with the claims in [59| and here we 
extended it to the case of self-gravitating objects, as long as we ignore finite size effects. Moreover 
we showed that we can adapt the transformation between NW and covariant SSCs found in 
for a spinning particle in a flat background to the (gravitationally bound) binary system. 



58| 



3C 



31 



331, 



51 



52] will be used to 



In a separate paper the results reported here and in [28l . 
compute the spin-orbit and spin-spin contributions to the energy flux and phase evolution to 3PN 
order [53 1. 
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Appendix A: Spin redefinition into precession form 



In this appendix we sketch the steps towards the precession equation. For the sake of notation, 
below it should be understood v — )• vi = vi/wj' = vi + . . ., and a — )• vi = aj + . . ., see Eqs. pTj 

MM- 

Moreover, all the equations will refer to the dynamics of particle one, even though V2 enters in 
the acceleration as in Eq. (J17p . In the spirit of abbreviation let us start by rewriting Eq. (j58p as 



A = a + vxa; + 6 A, 



(Al) 



where w = o;^^ = wq + ^v x ao at LO, with wq = uj^'~' (see Eq. [60]) ) and 



(5A = -v2ao + -v(vao). (A2) 



We will return to this last piece at the end. Let us start with the redefinition (same as in 



31 



33|, 



m 



S=(l-yjS + iv(vS), (A3) 

and move slowly. Taking time derivatives of this expression we have 

^ = _(v . a)S + [l - ^] S + ^a(v • S) + lv(a • S) + \{^r ■ S), (A4) 
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that we can write in terms of S after inverting the spin redefinition in Eq. ()A3 



S = — (v X S) X a (a X S) X V + tJo X f S v(S • v) ) + (v x S) x (a + v x a;) (A5) 



1 - - v2- -1 

— V V • [a;o X S + (v X S) X (a + V X i:;;)] — ((v x S) x a + (a x S) x v) + -(a x v) x v(S • v). 



Let us drop the bars from now on and consider only the relevant terms to NLO (recall S ~ v^S), 
then 

S = Wo X S H — (v X a) X S ( wq + -v x ap ) x v(S • v) — S • (v x a;)v (A6) 

1 v2 
H — V [(Jo X S + (v X S) X ao] • V ((v x S) x ao + (ao x S) x v) . 

Now we use [(v x S) x ao] • v = [(v x ao) x S] • v, and split it in two halves 



S = c^o X S H — (v X a) X S ( ijJq H — v x ao ) x v(S • v) — S • (v x a;)v 



1 
H — V 

2 



^0 + ^(v X ao) ) X S 



(A7) 



V V 

V - — ((v X S) X ao + (ao X S) X v) + - [(v x ao) x S] • v 



next we identify w = cjq + ^v x ao, and we obtain 



S = Wo X S H — (v X a) X S u x v(S • v) — S • (v x a;)v 



(A8) 



1 v^ V 

H — V (o; X S) • V [(v X S) X ao + (ao X S) X v] + - [(v x ao) x S] • v. (A9) 



Hence we realize 



1 v-^ V 

S = a;oxS + -(vx(a + vx uj)) x S — — [(v x S) x ao + (ao x S) x v] + — [(v x ao) x S] • v 



+ (v X S) X 5A, 



including the term in Eq. l\A2\i . In order to have spin evolution with constant norm, our remaining 
task is to handle the last four terms. Using Eq. (IA2p . these are 



v^ V 1 o 1 

6S = —— [(v X S) X ao + (ao X S) X v] + — [(v x ao) x S] • v H — v (v x S) x ao H — v • ao(v x S) x v 



v^ V 1 

— [(v X S) X ao - (ao X S) X v] + — [(v x ao) x S] • v + -v • ao(v x S) x v. 



(AlO) 
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which can be grouped as 

v^ V 1 

6S = — (v X ao) X S + - [(v X ao) X S] • V + -V • ao(v x S) x v. (All) 

Notice that the first term has aheady a precession form and, indeed, precisely completes the terms 
in the expression for A in Eq. (|Aip (the piece proportional to v in Eq. ()A2p cancels out in the 
cross product), so that 

1 V 1 

S = wo X S + -(v X A) X S + - [(v X ao) X S] • V + -V • ao(v x S) x v. (A12) 

In order to remove the left over terms we need to perform a 2PN shift. Garnering all the ingredients 
it is now straightforward to show that the following complete transformation (where we restore the 
bars, see Eqs. ([611 [751) ) 

removes the unwanted pieces in Eq. (IA12J1 and introduces one extra term given by Eq. (I75p . as in 
Eqs. dMllMD. 
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